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Abstract 

While the idea of gravity as an emergent phenomenon is an intriguing one, httle is known about concrete implementations 
that could lead to viable phenomenology, most of the obstructions being related to the intrinsic difficulties of formulating 
genuinely pregeometric theories. In this paper we present a preliminary discussion of the impact of critical behavior of certain 
microscopic models for gravity, based on group field theories, on the dynamics of the macroscopic regime. The continuum 
limit is examined in light of some scaling assumption, and the relevant consequences for low energy effective theories are 
discussed, the role of universality, the corrections to scaling, the emergence of gravitational theories and the nature of their 
thermodynamical behavior. 
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1 Introduction 



General relativity is often viewed as a low energy effective field theory [T], pretty much as a sort of "Fermi-like- 
theory" for geometrodynamics. This is largely due to its incompleteness at the classical level (due to singularity 
theorems), to power counting reasoning, suggesting that the Einstein-Hilbert action, with or without a cosniological 
constant, can be just the lowest order term of a series of terms, and, finally, to the fact that as a quantum field theory, 
the theory of an interacting spin-2 massless particle suffers from the same problems of four-fermions interaction^^. 

In the exploration of the physical content of general relativity, several clues have been collected suggesting that 
spacetime and its geometry might be emergent concepts in the semiclassical limit of an unknown theory, possibly 
pre-geometric in nature. On the other hand, all the quantum gravity models that we have at our disposao are based 
not only on rewriting general relativity in terms of more suitable variables, but typically imply a deep revision of 
the very concept of space and time [5]. 

It is often stated that the complete theory of quantum gravity will be able to cure the known diseases of 
perturbative quantization, as well as to provide an explicit computation (at least in principle) of the effective 
couplings determining the corrections to general relativity at increasingly high curvaturqj. Some results are known 
in this directions (let us mention, for instance, corrections to supergravity from string theory and effective equations 
in loop quantum cosmology [7]), but many other approaches, most notably spinfoams [8j HI [TOl [11] and group 
field theories (GFT) [121 ESI HI], are still rather far away from the stage in which concrete computations of 
physically relevant observables can be computed in a systematic way (see however, the attempts in |15[ I16j within 
a cosmological context). 

Many models have been proposed, but none has been proved decisive, to date, the reason of which can be traced 
back to the difficulty of extracting quantitites that could be cast in the form of physical transition amplitudes or 
cross sections, that can be concretely computed and compared with experimental observations. A rather remarkable 
exceptior0 comes from Causal Dynamical Triangulations (CDT) for which many numerical results of direct interest 
are available. We refer to [T71 [TSl for reviews, references and recent results. These last models, in particular, 
will be source of inspiration for many of the ideas that we are describing in the following. 

Parallely to these approaches, an increasing amount of work has been done in the direction of understanding 
gravitation as an emergent phenomenon. While the idea is not new, recently there has been a rather intense 
activity in the investigation of models based on condensed matter systems which display, in some regime, the 
emergence of effective metrics different from the one of the lab. On one side, the activity has gone in the direction 
of understanding better the properties of quantum fields in curved spacetime (Hawking effect, particle creation in 
cosmology and in strong fields, etc.) which is an important area that so far lacks experiments under controlled 
conditions. On the other side, there are some attempts to describe the properties of emergent dynamics for gravity, 
at least in very specific and limited situations. 

All the models that have been described so far do make reference to some background metric structure, thus 
making their impact on quantum gravity research rather limited. Indeed, as a consequence of the background 
structures, Lorentz invariance and diffeomorphism invariance are often broken symmetries. This has an impact 
both at the kinematical level, with a naturalness problem accompanying the appearance of Lorentz symmetry at 
low energy |20j . and at the dynamical level, with the graviton acquiring extra modes whose coupling to matter 
fields cannot be controlled in such a way to respect the observational constraints. To keep the idea alive, then, 
there are two possibilities. Either some specific mechanism is found to get rid of the problems encountered so 
far, or it is necessary to move forward and study completely different classes of models, where all the potentially 
dangerous features can be satisfactorily controlled, at least in principle. 

Recently, the group field theory (GFT) approach [T2j[13l[14] has raised a lot of attention as a possible realization 
of a (statistical mechanical) pregeometric theory that allows a background independent formulation of the path 
integral approach to quantum gravity. In this sense, we believe that this class of models do represent the ideal 

^On this, it is important to mention that, in recent years, this point has been reconsidered within the framework of asymptotic 
safety 

^With the exception of loop quantum gravity 4 , that takes very seriously the classical theory and proposes a Dirac quantization of 

it. 

■^However, transplanckian features [6] not necessarily related to planckian curvature regions might be difficult to grasp. 
^The status of string theory in this respect is rather controversial, given that the contact with effective field theories is more 
transparent, but still the correct limit including the standard model has not been identified. 
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arena for testing many of the ideas of emergent gravity |Sl[ll[23[Ml[13[13[271[Ml[inil3ni[SIll31[33l31]- One 

of the purposes of this paper is to elucidate the potentialities of such a framework in this perspective. 

One of the main point of GFT is that semiclassical continuum limit should be found in correspondence to a 
suitably defined critical point of such a theory. While many results have been accumulated so far ([351 133 133 
|38l[39l|40l|4ll|42l|43l|44l|45l|46l|^ just to mention some of them), still we lack explicit results concerning the 
presence of a critical point for the partition function in general models. However, in [15] the appearance of critical 
behavior for a class of models (that we will be interested in, incidentally) have been established, thus providing a 
strong basis for the discussion presented in the rest of this paper. 

In this paper we would like to address the general idea of finding imprint, on low energy physics, of the continuum 
limit defined as a suitable phase transition. Of course, lacking a definite calculation describing the general critical 
properties of GFT we have to rely onto some assumptions, that are somehow arbitrary but reasonable enough to 
describe a large class of possible (and plausible) critical behaviors. 

The aim we have in mind is actually twofold. On one hand we will try to pinpoint some very specific effect, 
ideally the pattern of the modifications to the low energy action for gravity. On the other hand, we will tackle the 
more general problem of the map between a microscopic pregeometric theory and a suitable effective field theory. 
While the results will be only partial, at this preliminary stage, they do shed some new light on the topic. Indeed, 
while we cannot yet give an answer to a physical question in terms of a number, we do propose a specific scheme 
to follow, in order to translate some formal results of the underlying model into a concrete effective field theory. 
The choice of GFT is by no means a reduction of generality: we believe that many of the comments that we will 
make can be exported to other approaches, even though the very structure of GFT allows to address some issues 
in a very direct way. 

Let us briefly sketch the scenario we are aiming at. The underlying pregeometric model (like GFT) will have 
the following properties: 

1. the description of matter and gravity is given in terms of pregeometric structures (matrices, tensorial objects 
etc.); 

2. the theory provides the partition function, as well as all the possible correlation functions of the theory; 

3. the theory provides the Schwinger-Dyson equations and the Ward identitie^ relating the correlation func- 
tions; 

4. certain composite operators will lead to the summatior0 over all manifolds, geometries and field configura- 
tions having some specified boundary data, corresponding to the various transition amplitudes we could be 
interested in; 

5. the microscopic theory depends on a number of free parameters, and the continuum limit is obtained by 
tuning the parameters to reach a critical point (hence a phase transition). 

In the language of critical phenomena, general relativity would be the result of a special phase transition where 
the order parameter is an entire three-geometry (i.e. a pair {^^,^ g)) (or an entire four-geometry, depending on 
which formulation of the theory one wants to describe), while the Einstein-Hilbert action can be regarded as a 
Landau-Ginzburg action for the order parameter (for a presentation of these ideas, see [49l [50l [51]). Of course, 
this description holds whenever we are well inside the geometrical phase, where the fluctuations in the metric are 
negligible. Near the phase transition, approaching the pregeometric phase, one has to take into account more 
details, and the mean field description of gravity must be abandoned. 

When trying to go for the statistical mechanical model underlying general relativity, one has to deal with 
the implementation of the symmetries and the structures characterizing the low energy effective theory (GR and 
standard model): diffeomorphism invariance, local Lorentz invariance, matter fields and internal gauge symmetries, 

^In the rest of the paper we will refer to the Schwinger-Dyson equations implicitly including the Ward identities and all the possible 
relations among the correlation functions that are implied by the properties and the symmetries of the partition function. 

^In general, we can expect that the sum is not really over the classical geometries etc., given that we know that this procedure 
would not really make sense. Rather, we might think that the sum is restricted to certain structures that somehow make the sum well 
defined (even just formally), for instance, simplicial complexes. 
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matter couplings and the equivalence principle. While some results |52[ [551 [Ml HH |5S] are already known in 
the literature, many questions remain open. For instance, one should clarify the emergence of diffeomorphism 
invariance in any dimension, the presence of a Lorentzian signature instead of a Riemannian one, the inclusion 
of appropriate matter fields and the recovery of effective field theories, the possibility of having anisotropic scaling 
between spatial and time directions (as in [56]), the clarification of the continum versus the semiclassical limit, the 
emergence of area laws for entropy associated to horizons and, obviously, unitarity of the quantum theory. 

Concretely, in this paper we will start from the very simple situation of GFT describing pure gravity. Despite 
rather unphysical, this case already presents a number of interesting conceptual difficulties that must be addressed. 
The goal is to investigate the consequence of the presence of critical points for the partition function of gravity (as 
it is defined in this way), and in particular how it is possible to relate the critical behavior to a low energy/long 
range effective field theory, by examining the shape of the wavefunction and guessing the shape of the corresponding 
Wheeler-DeWitt equation. 

The results that we will obtain are essentially three. First of all, it is clear that the continuum limit is not 
giving automatically the semiclassical limit. Rather, the effective Hamiltonian operator that we will define will 
still contain all the quantum corrections (coming both from second quantization of the gravitational field and from 
third quantization effects) . Second, we will draw a road map to pass from the knowledge of the critical behavior 
of the microscopic model to the formulation of a long range effective field theory. 

At this level, to really discuss the equation of motion, we will have to deal with the precise form of the Schwinger- 
Dyson equations of the theory, and with the way in which they operate in the large volume limit. However, we 
will elaborate some strategies to effectively deal with the problem without necessarily passing from the microscopic 
theory (which is however a necessary step that can be postponed a bit, but not avoided). 

Despite this, we can still establish an expected but still important result: the emergent theory, i.e. the coupling 
constants of the long range approximation will be determined by the critical behavior (in our case just by the 
critical exponent, but more complicated situations can be envisaged). Therefore, this leads to the conclusion 
that, close to criticality, the memory of the microscopic model and its parameters can be washed away, the only 
relevant feature being the universality class of the model. In other words, universality seems to suggest that the 
microscopic parameters of the model cannot be immediately identified with some macroscopic counterparts, as it is 
sometimes suggested, even when the construction of the microscopic model attempts to follow as closely as possible 
the macroscopic action that one hopes to obtain in the large scale limit. 

In addition, we will also estimate the impact of the corrections to scaling (corrections that will contain the 
information distinguishing the different microscopic models having the same critical behavior) when moving away 
from criticality and, finally, we will propose a concrete definition of entropy in terms of the microscopic models such 
that it might be compared with gravitational entropy as it is defined for semiclassical gravity (black hole entropy, 
De Sitter entropy, etc.). 

Contextually, we get a third outcome from this analysis. What is rather surprising is that, despite taking as 
a case study a particular model for GFT that would seem to be completely blind to information about curvature 
(in a precise sense that we will specify later), it still seems that, on large scales, the effective dynamics can be 
non-trivial, at least in the cases we have considered in this paper. This might be a coincidence, or, as it seems 
plausible, it is a manifestation of universality: given that we are looking at the system at criticality, the precise 
details of the microscopic action might be washed out, and the final long range action (in an ideal world, general 
relativity) might be just the leftover at the critical poini0. 

The plan of the paper is as follows. In section [2] we will introduce the setting, with the relevant definitions 
and the needed assumption. We will derive the general form of the thermodynamic potentials written in terms 
of the macroscopic extensive variables, like the total volume. After discussing some generalities of the problem of 
reconstructing the Wheeler-DeWitt equation in section |3l we will define, in section HI some generating functionals 
that will be naturally associated to the presence of macroscopic boundaries, allowing us to define weighted sums 
over boundary geometries. Furthermore, as a consequence of the scaling assumptions we will be able to argue the 
shape of the wavefunction in terms of the critical behavior of the system. Sections [5] and [6] will be devoted to 
the reconstruction of the Wheeler-DeWitt operator necessary to develop an effective theory for gravity, while in 
section [7] we will give an estimate of the impact on the results of the corrections to scaling. After discussing some 

^We will expand on this in the conclusion, after all the evidences have been collected. 
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conceptual implications of these results on topics like the origin of the gravitational entropy (section , we will 
summarize the results in section [5J where possible developments will be proposed. 



2 The setting 

In this paper we will take seriously the idea that spacetime is the result of a particular microscopic system being 
close to criticality. We will go further with our assumptions and make an additional stipulation on some properties 
of criticality. The hypothesis we make is rather natural in light of what we know about critical phenomena. In 
particular, we will assume that the partition function for quantum gravity, thus defined, will display characteristic 
scaling behavior with respect to the parameters controlling the approach to the critical point. Of course, this 
should be accurately motivated only from the analysis of the underlying microscopic statistical mechanical model, 
but still this represents a natural option, in light of the body of work related to matrix models j57) . The recent 
result [IH] strongly supports this point of view. However, the discussion presented in the rest of the paper could 
be adapted (at least in principle), to any other form of critical behavior. 



2.1 GFT: actions and notations 

Let us take, to make the discussion concrete, the case of the partition function of matrix models/GFTs. This gives 
(even if only formally) the sum over all the possible simplicial complexes (with appropriate topological restrictions), 
each one weighted by an amplitude which is uniquely determined by the original matrix model/GFT action. 

A group field theory [T^l [131 HI] can be seen as a generalization of matrix models ^57^ and models like the ones 
proposed by Boulatov [SS] and Ooguri [SH] to the definition of the partition function of gravity in dimension three 
and four, respectively. The starting point is a field theory (for instance, complex scalar field theory) over a certain 
number of copies of the group manifold G, 

(1) 

where d is the number of the dimensions of the spacetime that we want to discretize. In the remainder, to avoid 
clumsy mathematical formulae, we will use the notation 

(t>a...b = <P{ha,--,hb). (2) 

The theory is defined by means of an action 

SGFT[(t)] = J {dhf(j)a,...aAai...a^ ^ Z / ^'^d)^ ^^^'^^'>'f'<^l--<---'^ai+\..a''/^ + ^C- (3) 

where V is a certain function of the group elements encoding the combinatorics of the d— dimensional symplex. 
Of course, more complicated choices are possible, enabling the implementation of specific properties of the class 
of histories that one is summing over. These important modifications will not influence the general gist of the 
discussion, and we will ignore them. The partition function of the model is defined as usual by: 

Z{g)^ J V^eM-SoFTm. (4) 

It is possible to work also with the partition function defined with exp(— iS'g^ft). However, this would not add 
much to the discussion. The very definition of such a partition function is far from obvious, and a suitable 
regularization procedure, defined in terms of a expansion |35l 1361 157] . must be adopted to provide a workable 
partition function. Very briefly, the partition function will depend on the coupling constant, g, used to deflne the 
perturbative expansion of the theory. The Feynman graphs associated to this perturbative expansion are can be 
described by simplicial complexes, while the amplitude of each simplicial complex appearing in the sum is nothing 
else than a spinfoam amplitude. It is easy to see that, taking 

W{g) ^ log Z{g), (5) 
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to eliminate the summation over disconnected graphs, we can define a partition function for d— dimensional quantum 
gravitjU. We refer to the literature for more complete expositions. 

The function W can be given in terms of its pcrturbative expansion, 



(6) 



r 



where the summation goes over all the simplicial complexes F compatible with the Feynman rules of the model, 
n(F) is the number of simplices of the complex and A{T) is an amplitude that depends on the complex as well as 
on its colorings, when relevant, that can be understood as a Boltzmann weight factor in the statistical sum. 

The function W can be used to compute certain expectation values. For instance, it is easy to realize that, 
when using its perturbative expansion (|6l), the expectation value of the number of simplices reads 



This quantity might be related, for instance, to an average volume, once the mapping between the combinatorial 
structures of the graphs and the metric information is established. Despite this is not obvious, we will do so, in 
the rest of the paper, as we shall specify in a moment. 

Symmetry requirements on (/) and on V can lead to further specifications about the class of simplicial complexes 
that are entering the sum, as well as their colorings with group representations and intertwiners which are encoding 
the geometrical properties (metric, connection, curvature) at the discrete level. Here we take a departure from the 
mainstream approach: while we take V{{h}) to be the correct one to lead to simplicial complexes of the appropriate 
dimensions, we do not use any projection on the field on the group manifold to enforce any form of invariance under 
special subgroups of G x ... x G. As a consequence, it will be impossible to store any form of geometrical data on 
the simplicial complex, beyond the combinatorial structure, as it happens in more refined models implementing 
spinfoam amplitudes. The only thing we will be able to count is the number of simplices we are using. 

With respect to other GFTs for quantum gravity, like the one of Boulatov, which imposes a certain symmetry 
on the field (f>, necessary to reproduce the amplitudes of the Ponzano-Regge model for three dimensional gravity, 
we might think that we are indeed working with a theory of triples of edges that are not necessarily closing into 
triangles. However, we could also look at the model as a model for equilateral simplicial complexes. Besides 
the information on curvature encoded into the deficit angles (ultimately related to the combinatorial structure of 
the model), there is no other possibility to weight differently complexes with different curvature structure. This 
represents a great limitation imposed on our analysis, which will seriously affect our possibilities to make serious 
contact with models for quantum gravity proposed so far. We will come back on this point later, showing that this 
might be, in fact, a welcome property. 

In analogy with matrix models, the continuum limit is thought to be achieved by tuning the coupling constant 
to a critical value, gc, to be determined by further investigations. The reason for this is that, at a critical point, 
the number (n) diverges, signalling that the "average triangulation" will contain an infinite number of simplices, 
a necessary condition to speak about a continuum limit. This is the reason why the understanding of the critical 
behavior is of paramount importance to understand completely the continuum, semiclassical limit of quantum 
gravity (in this particular approach). 

For the rest of the paper we will assume to have found this critical point. Close to a phase transition, the free 
energy (as a function of the coupling constant) can be separated into a regular, analytic part and into a singular 
part. All the critical behavior is controlled by the singular part, while the regular part contains the information 
about the behavior of the system away from criticality. 



Let us start from the assumption that at criticality, the logarithm of the partition function for GFT has a singular 
behavior of the form: 




(7) 



2.2 Criticality 



W = a{g - g,r (1 + .f{g)) 



(8) 



For a thorough discussion about the Feynman's path integral approach to quantum gravity, see |60| 
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where f(g) is a function of the couphng constant, analytic in a neighborhood of gc- Of course, this is a rather 
arbitrary assumption, to date, that consists in exporting verbatim the lessons of matrix models, at least in their 
simplest incarnations [57 . Of course, one can imagine different forms for the non-analytic part of W: only an 
explicit calculation can settle the matter. We will come back to the point in the concluding remarks. However, as 
said, the results of 48 support this conjecture. 

This stipulation ([S]) , which is the basic assumption for this section, amounts to a specification of the universality 
class of the critical point, determined essentially by the exponent 7. We assume that other contributions to W do 
not influence the critical behavior, while they represent the nonuniversal corrections away from criticality. These 
are precisely the kind of corrections that could allow us to discriminate between two distinct microscopic models 
having the same critical behavior. We will neglect them, for now, and therefore we will be describing only the 
general signatures of critical behavior, not the ones specific of the given model. 

Generically, it is expected that the exponent 7 could depend on dimensionality of the model, the kind of 
interactions, presence of additional degrees of freedom, pretty much in the same way in which it has been discussed 
in the context of matrix models. Define, for convenience, 

C = log 5- (9) 

If we want to estimate the average volume (or surface, if we consider matrix model), we need to compute the 
average number of vertice|f|, which is seen to be: 

, ^r ^ 1 dW 1 dW 

(^^)(9)=n=-g^^-^, (10) 

This number is the average number of elementary cells in the Feynman expansion. Clearly, g ^ gc corresponds to 
the limit of infinite combinatorial volume^. Notice that with one parameter at our disposal we cannot really probe 
more than this, at this stage. As usual, to decide whether this number is physical or not we need to see whether 
this is somehow related to the physical volume probed by matter fields. To uncover other properties of the average 
configuration, like topology, curvature etc., we would need more external parameters/coupling constants. 
In our case, as a consequence of (|8]) 

n= ^ ^, (11) 

where we are dropping the angular brackets to avoid clutter. We can invert these relations, to get: 

9= , C = log (12 

n — 7 \n — J J 

Notice that these relations holds only close to criticality, when 1. 

This number is just combinatorial in nature, and does not immediately lead to a physical volume, for which the 
metric structure needs to be established. To do this we need to introduce a length scale in the game, e and define: 

V = e'^n. (13) 

The nature of e requires further discussion. If we consider just g ^ gc we just get an infinite n limit, not a 
continuum limit. As it is well known, in the case of lattice field theory, taking an infinite number of lattice points 
will not lead to any continuum limit of the theory. To speak of the continuum limit we need to check that the 
correlation length of the field theory becomes macroscopic with respect to the lattice spacing. Therefore, to discuss 



Again, remember that in the perturbativc expansion, each vertex is dual to a d— simplex. Furthermore, the function W can be 
expanded into 

W{g) = J2 9"^^^Mn 

F; connected 

where n(r) denotes the number of simplices contained into the complex F, and j4(F) is an amplitude determined by the combinatorial 
structure of the Feynman graph (including symmetry factors). 

^"Sometimes it is mentioned that the GFT coupling constant is related to the cosmological constant. Indeed, as these expressions 
show, it is in fact related to a variable that is the thermodynamical dual of the average volume. 
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in a physically relevant way the continuum limit we need to introduce matter fields, representing rods and clocks 
which will give the precise mapping between the combinatorial quantities defined at the level of the Feynman 
expansion and the dimensional quantities of physical interest. In the setting we are working, this issue is left aside 
and it is implicitly assumed that all this procedure of introducing matter fields etc. is directly included in the 
appearance of the parameter e which translates the combinatorial volume (counting of elementary volumes) into 
a physical volume. We could imagine, for the sake of the discussion, that e is the ratio between the correlation 
length ^ (defined in combinatorial terms), and the lattice spacing a: 

e=|. (14) 

This ratio will represent the ratio between the Planck scale and the ruler constructed with the correlation length 
of the given field. With this in mind, we will drop the e from our expressions. 

2.3 Consequences 

Define the free energy to be: 

F^~\ogW (15) 
so that W is maximized when F is minimized. In our case 

F = —7log(.g — .gc) + const. (16) 

This is a thermodynamic potential of the system when ^ is used as an intensive macroscopic variable to describe its 
properties. On the other hand, we could define another thermodynamic potential, say fl, defined as the Legendre 
transfornf"] of F, where now the system is described by the conjugate macroscopic variable, n: 

n{n) = F(C(n)) - C(n)n = -7 log f - n log ( (17) 

Let us manipulate this expression 



n{n) ^ -7 log(7.gc) + 7 log(7i - 7) - n log(.gc) - n log (18) 

\n~-fj 

Now, as we said, we are interested into the limit in which n 3> 1, for which: 

Ti 

n — J n, ~ 1 (19) 

n — 7 

Here it is clear how 7 represents a sort of elementary volume. One would be tempted to use it to generate and 
expansion of in powers of 7. However for these corrections to be really important, we need n « 7, which is a 
regime we are excluding from our assumption of being near criticality. If we were considering a system which is 
too far away from criticality, we would need to include all the nonuniversal contribution to the partition function 
that we have neglected, here, from the very beginning, which presumably will give relevant contributions. We will 
describe them in section [T] For the moment, we will neglect all the terms of order 1/n. In this limit: 

r!(n) «7log(n)-nlog(5e) (20) 

where we have neglected a constant which does not depend on the thermodynamic variable n. 
We can use fl(n) to define a probability measure for the volume n, 

P{n) = — exp(— f2(n)) = — n'' exp(— an), (21) 



^^With the conventions commonly used in thermodynamics to define the procedure. 
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where 

a = -loggc, (22) 

and 7 is the critical exponent. 

The normahzation function S is easily computed 



/>oo 

"E. — I n'' exp(— an), 



(23) 



where niR is an infrared regulator needed to get rid of the small volume contribution, that is not within the domain 
of applicability for this model. Keeping this in mind, we will consider the case niR = 0, for convenience in the 
calculation but also because, for the case 7 > 0, the term automatically suppresses the contribution coming 
from small volumes. 
Therefore 

_^,u-< du r(7 + i) 



rC em{-an)dn ^ I due = —-!——!- (24) 

With this statistical measure we can compute average value for various quantities, like the average value of 
some power of the volume 

(„') 4 exp(-„„M„.-lI<l±|±ll (25, 

This expression can be further simplified using the recursive properties of the Euler's gamma function 

fc 

V{z + l) = zT{z), V{z + k)^ \\\{z + i)\T{z) (26) 

to give 

(«') = «-MTT(7 + 0) (27) 





Define the average volume to be 



Clearly, for consistency 



T + 1 

r = (n) = ^— . (28) 
a 



> 1 ^ a < 1 (29) 



which means that gc must be very close to one, so that the most probably a weak coupling expansion cannot be 
used at all. This is consistent with the fact that, to have arbitrarily large Feynman diagrams we cannot truncate 
the perturbative expansion at a given finite order. 

Similarly, we can define the spread of the distribution, by considering 

(ri^) - y2 1 

^ = I., = (30) 

Equations (j29p and (j30[) suggest that we can replace a and 7, which are defined by means of the microscopic theory, 
with some corresponding macroscopic quantities which control the mean value and the typical fiuctuations around 
the mean value. In particular: 

7=^-1, a = i (31) 

whence 

r(A-i) / 1 n 



P{n) = ^-^n- -^exp(^---j (32) 

The way in which A enters this expression suggest to think about it as a "temperature" . Nonetheless, this 
identification is only formal. Interestingly, A is determined by the critical exponent 7 which describes the way in 
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which the continuum Hmit is approached, and hence cannot be varied ad hbitum without varying the (universahty 
class of the) underlying microscopic model. Introducing the notation 




(33) 



our expression for P becomes: 




(34) 



Notice that the requirement that the resulting distribution has a very small spread around the mean value (which 
in light of our interpretation of this kind of models corresponds to a spacetime which is essentially classical) , A <C 1 , 
which implies that 7 3> 1. In terms of the "inverse temperature" /?, this condition can be seen as the requirement 
that the ensemble is kept at a very small temperature (but again, this interpretation is only formal, since the critical 
exponent 7 is not a parameter which can be easily controlled, being determined by the underlying dynamics of the 
model) . 

3 Wavefunction and the reconstruction of Wheeler— DeWitt 

The discussion of the previous section was focused on a very primitive approach. While for models of two dimen- 
sional quantum gravity the use of equilateral triangulations proves to be rich enough to get some nontrivial results 
(see, for instance |57j). for higher dimensional structures we have the additional problem of introducing, in one 
way or the other, a notion of curvature at the microscopic level. For instance, in GFT the gravitational action 
can be introduced by certain projections on the group manifold that one is considering, in such a way that, when 
developing the perturbative expansion of the theory, each vertex of the Feynman graphs is associated to a suitable 
amplitude which encodes, in the appropriate language, a discretized version of the Einstein-Hilbert action. For 
instance, in the case of three dimensions, the vertex is associated to a Wigner's 6j symbol, which, as proved by 
Ponzano and Regge [HI] is deeply related to the cosine of the action used in the Regge calculus approach to classical 
(and quantum) gravity 

In the case of causal dynamical triangulations for four dimensional quantum gravity, besides the volume term 
counting the number of 4-simplices, the curvature is introduced by appropriately weighting the 2-simplices on which 
curvature is localized (in a distributional) sense. There, one has two distinct parameters, one for the curvature and 
one for the volume, to control the system, to drive it to criticality and to obtain different phases. We refer again 
to the literature [T71 [TH] for more details and references. 

In the case of models like GFT, the only parameter at our disposal is the GFT coupling constant, with the 
interaction term encoding the entire gravitational action. Therefore, it seems difficult to disentangle the various 
contributions from curvature and cosmological constant already at the level of the partition function. Therefore 
we need something different. An option is to consider certain correlation functions of the GFT model which 
might have a more direct physical interpretation for quantum gravity. In this perspective, there are two natural 
candidates: the "wavefunction" for quantum geometrodynamics and the transition amplitudes between different 
boundary geometries. 

With the term wavefunction we will denote the particular correlation functions of the GFT model corresponding 
to the summation of triangulations possessing a single boundary with given data, representing a discrete geometry. 
The transition amplitudes are those particular correlation functions that correspond to the summation over trian- 
gulations possessing two boundaries with (in principle) two different sets of data encoding geometrical properties. 
Examples of transition amplitudes have been considered in the past, especially within the realm of the spinfoam 
approach to quantum gravity, and lately to quantum cosmology [15'. On the other hand, the wavefunction ap- 
proach originated already with DeWitt 63 , in an attempt to define canonical quantization of gravity based on 
the implementation of the constraints of general relativity as functional differential operator acting on a function 
defined over the configuration space of canonical general relativity, superspace. We will follow this second point of 
view, in particular in light of the Hartle-Hawking proposal |64) that fits pretty naturally our framework. 

The Hartle-Hawking proposal is a prescription to define the wavefunction of the universe in terms of a Euclidean 
path integral defined over compact geometries having a unique boundary whose intrinsic geometry is the particular 



9 



point of the superspace where the wavefunction is evaluated. Clearly, our choice to work with the partition function 
defined with exp(— 5) instead of exp(i5) fits particularly well this scheme. 

There are interpretational issues on what we should define as wavefunction. In ordinary quantum mechanics, 
a wavefunction is usually defined as a projection of the quantum state on an orthonormal basis constructed out of 
some complete set of commuting Hermitian operators, associated to certain functions of the classical phase space 
of the theory. In the case of gravity the situation is more complicated, given that the natural domain of definition 
of a wavefunction, superspace, needs to be put in correspondence with operators whose eigenstates are, indeed, 
3-geometries. Even if we would succeed in providing such a construction, there would be the difficulty of relating 
these states with the ones that we could typically construct by means of discrete structures. In this respect, the 
most refined and rigorous results are obtained within the formalism of spin networks. 

Again, let us take the case of GFT and the wavefunctions that we can construct out of them. Let us focus on 
the four dimensional case. The quantities that we can construct (and perhaps compute within some approxima- 
tion scheme) are some correlation functions corresponding to three dimensional simplicial complexes, with given 
topological properties, colored with some representation label, say j if we consider SU{2). Let us take the point 
of view that these quantities are really possible states of quantum gravity, solving some discretized version of the 
Wheeler-DeWitt equation. Let us denote them as |A,{j}). Clearly, what we would need to define is a scalar 



which would allow to make a precise connection between the discreet and the continuum formalism. This is not 
just a side issue. Take the case of a single line of length L: we could discretize it with ni segments of length 
L/rii or with n2 segments of length L/n2. Both these discretizations are legitimate, a priori, and hence they might 
contribute, with different weight, to the definition of the state that we call "a macroscopic loop with length L". 
The problem of the coarse graining of microscopic geometrical information encoded into graphs, and the definition 
of states which have definite global geometrical properties is still subject of active research |65) . 

Furthermore, as we have already said, to be really able to address physical questions, one has to introduce 
matter fields, that can be used to introduce notions like physical time (as it happens, for instance, in the case of 
quantum cosmology [66]). The realization of group field theories containing realistic matter fields in their Feynman 
expansion is still work in progress (see, however, |llj). In addition, we have to include the notion of continuum 
limit/criticality, and, finally, of semiclassical limit, which will involve the definition of suitable quantities and 
parameters associated to these procedures. 

While this constitutes a whole research program that cannot be exhausted in one paper, we can still try to make 
conjecture and work out some of the consequences of them, pretty much within the lines of what we have discussed 
in the previous section. In particular, we will implement, in a different way, the scaling hypothesis that we have 
made in the previous section, now treating the three geometry as a thermodynamical variable. Let us be more 
specific on this. While the complete specification of a geometry might be a hopeless task, we might imagine that we 
can treat situations with high symmetry so that many degrees of freedom can be effectively neglected. Therefore, 
we assume to restrict the analysis to describe situations in which we can completely identify the geometry that we 
are interested in by means of a certain number of functionals, for instance the total volume, the mean curvature 
etc. For the sake of the argument let us assume that we can restrict this infinite family of functionals (local and 
nonlocal) to a few of them, say the total volume and the average scalar and square curvature 



where we are fixing once for all a background manifold and we are restricting q to be of a certain family, such 
that these three coordinates give a complete specification of which element of the family of metrics q we are talking 
about. This approach has all the limitations of the symmetry reduced systems, but still will give us the possibility 
of concretely draw a road map. 

Our wavefunction will be constructed out of these quantities 




(35) 




(36) 



(37) 



Or a suitable projection. 
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and, in principle, its shape can be computed on the basis of suitable computations made at the microscopic level 
[67] . For now, we will assume that this has been done, and that the wavefunction (j37p is known explicitly as a 
function. By the very construction, this quantity obey the diffeomorphism constraint equations: 

V. (^^m) = 0, (38) 

given that ip is constructed with diffeomorphism invariant quantities. The true question is: what is the Hamiltonian 
constraint such that this wavefunction solves the corresponding Wheeler-DeWitt equation? As a subordinated 
question, what is the semiclassical limit of the theory obtained in this way, and finally, what are the quantum 
corrections that we can expect, with their associated phenomenology? Such a question is nontrivial even in more 
traditional approaches to the path integral for quantum gravity |68i I69| . 

The situation is clear: a formalism like GFT, or spinfoams, promises to compute quantities like or more 
general matrix elements. However, to really make contact with the language of effective field theories, which is the 
language in which we can talk about experimental tests of theories, these days, we need to translate this quantity 
in terms of Lagrangian, or, equivalently, differential equations. 

As we said, the theory that allows us to compute the wavefunction should also provide some relations between 
the wavefunction computed at different points. This is the point at which the Schwinger-Dyson equations for 
the correlation functions come into the game. These equations will be, generically, recursion relations between 
correlation functions with different number of points and with different colorings. When performing the continuum 
limit, one has not only to compute the continuum limit of the correlation functions, but also the continuum limit of 
the equations relating them. This is what has been done, within the context of matrix models, with loop equations 
[571 1701 171) . Loop equations are recursive equations relating matrix correlation functions associated to loops with 
different lengths, and, ideally, in the continuum limit, they should give rise to the equations for the continuum 
dynamics [71 [71 [71]. 

For the case at hand, the situation is still not clear. It is plausible that, given the relation between the 
constraints of GR and diffcomorphisms, one should use the Ward identities associated to the symmetries realizing, 
at the microscopic level, diffcomorphisms. To this purpose, the work in [44] is of special relevance. Nonetheless, 
still much more is needed, especially in higher dimensions. 

Without the Schwinger-Dyson equations, or some recurrence relations among correlation functions, the formu- 
lation of an effective field theory becomes rather difficult and subject to uncontrolled arbitrariness. In this paper 
we will able to generically address this arbitrariness, but a solution can be obtained only when the exact statistical 
properties of GFT will be elucidated. 

In the remainder, we will try to extract as much as possible from the only knowledge of the critical behavior, 
without having at our disposal any information about the dynamics in the critical limit. While this represents a 
limitation, we will still be able to uncover a number of interesting properties. 



4 Generating functions 

A key tool to define consistently the wavefunction, within a combinatorial approach, is by means of a suitable 
generating function. This will be done by adapting the procedure used for matrix models to the case of higher 
dimensional spaces. In the case of (single) matrix models, it is customary to define the generating functions 

where (...)c denotes the expectation value obtained counting only the connected diagrams, x,y are real numbers 
and M is the matrix for which the statistical system is defined. These particular correlation functions represents 
summation over two dimensional triangulations with single or multiple boundaries, given by loops whose average 
length is controlled by the "fugacities" x, y, representing the thermodynamically conjugated variables to the length 
of the loops. In other words. 
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where £ is the length of the loop. This result might be derived in the same way in which we have counted the average 
number of simplices from the derivative of the partition function W . Furthermore, these correlation functions, as 
a consequence of Schwinger-Dyson equations, obey a hierarchy of equations |57] . known as loop equations, which 
represent a direct path to the understanding of the continuum limit of the models [71| . 

The reason why this approach should be preferred is that it explicitly addresses the question about what super- 
positions of different microscopic configurations (i.e. boundary lengths) will contribute to the state macroscopically 
labeled by a certain geometrical property {i.e. the macroscopic length). Clearly, different choices of generating func- 
tions will generate sums that will weigh in different ways the various contributions, representing different statistical 
ensembles, according to the number of parameters introduced. 

In this section we propose a generalization of these generating functions to the case of higher dimensional 
models. The objective is to define in a consistent and controlled way the sum over different discretizations of the 
same boundary. For our purpose, we will need to consider just the case of a single boundary, but we will envisage a 
way to define the generating function for multiple boundaries. Let us start with the general structure of the result 
that we want to obtain. We want to define in a consistent way a weighted sum of correlation functions which have 
specific combinatorics associated to a specific combinatoric structure for a certain D-1 manifold, something of the 
form 

G(A) = J2 A"(^°-^)G(r,3_i) (41) 

Td-i 

where A is a parameter needed to control the average size of the boundary (i.e. a chemical potential dual its volume) 
and G(r) a GFT correlation function associated to the boundary triangulation F^i-i. 

The crucial difficulty in defining a sum over all the possible discretizations is to generate them. However, for 
this purpose we could use the fact that with an auxiliary GFT associated to a — 1 dimensional model we can 
accomplish this tasll^. Let us work, for instance, in three dimensions. Then, define 



G{X,9) ^ J 'Di(;T>(l)exp - Q J {dgfipabAa + ^ J {dgfilJabil^bci^ca(f>hih2h3 ] + 



'^[{d9)'M.ba+'^lMaefMbf 



(42) 



where hi = gaQb, ^2 = 5h5c, ^3 = <7c5a- The first term in the exponential generates a perturbative in A which would 
consists of all the possible two dimensional triangulations, if instead of the 4>abc we would have just 1. The presence 
of ci 

abc implies that the Feynman expansion will be in fact a sum of correlation functions for the three dimensional 
model, i.e. a sum of amplitudes of random 3D triangulations with random boundaries. 

By modifying the properties of the auxiliary group field theory (presence of colors, different combinatorics, 
modified kinetic and potential terms), we can control the class of boundaries we can sum over. The study of this 
generating function goes beyond the scope of the paper. However, it represents a new direction for the study of 
these kind of models of quantum gravity and deserves further study both for its definition and for the investigation 
of its properties. 

Clearly, the definition given in three dimensions can be exported to the case of higher dimensions. Also, by 
inserting other auxiliary group field theories, we can define generating functions which are associated to summation 
over simplicial complexes with multiple boundaries. 

The definition (j42p provides a summation over three dimensional simplicial complexes with arbitrary two di- 
mensional boundaries, and the parameter A represents the thermodynamical dual to the size of the boundary (the 
only variable that at the moment we are able to control). 

The usefulness of this generating function is clear: it allows us to define consistently the sum over all boundary 
triangulations, once we specify the statistical property of the sum, controlled by A and g. In particular, we can 
conjecture the existence of a critical behavior that corresponds to a continuum limit in which both the boundary 
and the bulk are triangulated with an infinite number of blocks (the caveat on the definition of the continuum limit 
in terms of appropriate correlation lengths, discussed in section II, still applies). As such, we can imagine that, 
near the criticality, the connected part of G, 

Q(A,5) = -logG(A,5) (43) 



'An analogous idea, applied to a rather different aspect of GFT has first been suggested by F. Hellman |75| . 



12 



could have a behavior similar to the one of the free energy W studied in section II, namely 

Q(A, g) « (A - XMfia - 9c{X))\ (44) 

were we have included the possibility that the critical values Ac, 5c (and, in principle, also the critical exponents) 
have a dependence on the value of the other coupling constant. For the moment, we will drop this possibility, but 
we will come back to it in section [T] 

As a consequence, we can proceed in perfect analogy with section[2]and get to the same results, once we replace 
g with A. In particular, we can infer that a boundary with size A will have a probability 

aS poo 

P(A) = — exp(-a'A), e= / A^eM-a')AdA (45) 
" Jo 

where 

a' = - log Ac (46) 

Clearly, this formula is valid only in the continuum limit and in the case in which we are very close to the critical 
point (Ac, gc), and neglects a lot of possibilities for the critical behavior. In this respect, a direct calculation cannot 
be circumvented. However, given that to date we do not have it, we will try to get as much as we can from this 
conjecture. In the rest of the paper we will try to extract as many consequences we can from equation (|45p. We 
will come back later on the issue of corrections to scaling. 



5 Reconstruction of WDW 

Now that we have collected some ideas concerning what kind of behavior we should expect, let us come back 
to the problem of translating a wavefunctional in superspace in terms of a Wheeler-DeWitt equation. One of 
the interesting question raised by the previous analysis is: what is the effective (field) theory underlying the 
results obtained via GFT? In the investigation of (causal) dynamical triangulation, the issue of the formulation 
of an effective theory describing the various phases of the model is of great importance. There, the presence of a 
microscopic partition function that allows a more or less direct connection with a transfer matrix or an Hamiltonian 
allows us to go as far as the construction of an effective fiducial Hamiltonian, at least in certain simplified cases 

m- 

We assume that, due to the nature of spacetime as a continuum limit, the microscopic system is at criticality, 
and that this fact manifests itself into the scaling properties of the correlation functions, in terms of the macroscopic 
variables: 

V'M = (/iM)"exp(-/3(/i(g))") (47) 

as an example. Assume, for the moment, that we further reduce our superspace to the space of spheres, and that 
the only free parameter determining totally the geometry is the radius of the sphere, r. Then: 

V'[q] = (r)3^exp(^-|r3"^ (48) 

where we have changed slightly the notation to include the information that what does matter is the action of the 
system, weighted in unitqlj of H. Furthermore, we have introduced 3a to keep in mind that, in 3 + 1 dimensions, 
when a — 1 we recover the standard large volume behavior of the wavefunction. We can now try to make contact 
with what we know. We want to find a linear differential operator L, such that 

L^Pir) = 0. (49) 

Moreover, we know the semiclassical limit, the equation, treated in the WKB approximation, should lead to an 
eikonal equation which is nothing else than the Hamilton- J acobi equation for the classical limit, that is general 

^^The emergence of units is indeed a big foundational problem that needs to be addressed. 
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relativity. As a consequence, we can cook up a linear (self-adjoint) differential operator with at most two partial 
derivatives with respect to r 

ixriai = h-^,+^hu{r)^^ + (v{r) - '-^{r)) (50) 



Qj.2 



dr 



2 dr 



and with a number of free functions to be determined. 
In this case 



whence: 



Qj.2 



37(37-1) ^ 9«^/3%6c.- 

y,2 ^2 



2 _ ^j,3a-2 



(51) 



(52) 



Now, if the wavefunction is real (as it happens for ground state functions of one dimensional systems), and in this 
parametrization this means /3 G M, the function u[r) is identically zero while the function v(r) is simply 



v(r) 



2 o2„6a-2 



3a-2 



37(37-1) ,2 



(53) 



Obviously, the theory will not be truncated in an expansion in /i, given that the starting point is not a semiclassical 
limit, but rather a large volume continuum limit in which all the quantum corrections are sistematically included. 
In this sense, this shows how the critical behavior encodes just the continuum limit, and not the semiclassical one, 
which should be achieved otherwise. Let us stress again that h in itself needs to be defined, as the reference action, 
and that the models, being purely combinatorial, need to be manipulated in order to identify, at criticality, the 
emergence of scales and physical units. For a similar case, involving the emergence of the speed of light, see ^5]. 

The operator found in this way can be used easily to develop a systematic effective field theory and to make 
definite predictions about physical observables, rather than using the wavefunction itself^. In principle, knowing 
the effective potential in the WDW equation we might want to go backwards and define a classical effective action 
leading to it in the suitable limit (perhaps involving some effective order reduction, as in |77[ 178] ). 

Extrapolating the discussion from the mere Riemannian quantum gravity models, it is clear that in the way in 
which we have proceeded, scaling is limited to the dependence on spatial properties, and the theory emerging needs 
not to be the 3+1 version of a Lorentz invariant theory (especially when matter fields are included). Consequently 
there must be a key argument to keep Lorentz symmetry under control. This should be provided by the Schwinger- 
Dyson equations and all the other relations among the correlation functions, which are enforcing at each level the 
microscopic symmetries (for instance, diffeomorphism invariance of the macroscopic theory as encoded in terms of 
a symmetry of the microscopic GFT). 

The fact that Lorentz invariance must be implemented already at the microscopic level seems to be rather 
compelling in light of the evidences that have been collected in the latest years about effective field theories with 
Lorentz violation effects systematically included [M]. However, the realization of the symmetry at large scales 
is by no means guaranteed, since there might be some form of spontaneous symmetry breaking which cannot be 
addressed without the complete understanding of the continuum limit. 



6 Characteristic curve method 

There is another way to get a better grasp of the inverse problem of finding the Wheeler-DeWitt equation given 
a solution, and it uses the method developed to solve first order partial differential equations, the characteristic 
curve method. The key tool we would like to exploit is the relation between the theory of linear partial differential 
equation, the WKB method and the Hamilton-Jacobi formulation of classical mechanics. 

The idea is simple. If the dynamical equation of motion has the form of a differential equation (of arbitrary 
order), it is possible to find approximate solutions by means of the cikonal approximation, which consists of 
approximating the solution with a function of the form exp(iS'(x)/L), where L is a small parameter (in appropriate 

^^On this point, see the comment in the introduction of chapter 4 of |79| . 
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units). While the method is used extensively in quantum mechanics, specialized to the case of second order PDE, 
its validity is more general [80 . 

Within this framework, given the Hamilton function, we can find the Hamiton-Jacobi equation it satisfies, trying 
to go backwards to uncover the nature of the differential equation which is ultimately its origin. The method, as 
we will see, will not be free of ambiguities, but still they will be drastically reduced, with respect to the very hand 
waving analysis of the previous section. 

As it is well known [ST , to solve a first order PDE of the form 

H{dS, x)=0 (54) 
where x E U C M", 5 : C/ — > R, one can proceed by defining 

P^ = d,S (55) 

and hence the canonical system 

. dH . dH ■ . 

With this method it is possible to generate solutions of the PDE obeying the desired initial conditions. Now, we 
are interested into the problem of finding H, once S is given. For instance, 

S{x) = fix) (57) 

with a certain, given function f{x). Knowing S we can compute p{x) — f'{x) and hence the task is reduced to the 
problem of finding a function H such that 

H{f'{x),x)^0 (58) 

Clearly, there will be more than one function satisfying this condition, and we will be left with this fundamental 
ambiguity that needs to be solved with other means. 

Instead of giving some abstract general derivation, we will specialize to the cases in which we are interested, 
that is to say, to functions of the form 

f{x) — log(a;''' exp(a;'')) — x'' ~ 7 log a; (59) 
which is the sort of functions we could be interested in, with x being related to some form of average size. Clearly: 

p(x) = bx^-^ - ^ (60) 

X 

On the other hand: 

whence we obtain a PDE for the function H{p, x): 

dp OH 9H _^ 
dx dp dx 

or, using the expression for p: 

This PDE can be solved via the characteristic curve method. In particular, we can define a new canonical system: 

It is straightforward to realize that the function H is constant along the curves: 

p - bx^~^ + - = const. (65) 

X 



15 



which means that: 

H{x,p) = $ - bx^-^ + ^) , (66) 

with $ an arbitrary function. To fix we still have to use the information coming from the SD equations, a step 
that cannot be completely avoided. However, the ambiguity in the system is entirely contained in the unknown 
function and hence, by making appropriate assumptions, one could develop a systematic effective theory for the 
classical limit. 

We can exploit this technique, then, to find the classical Hamiltonian which is associated to the WKB approx- 
imation of the solution of the equation we do not need to solve (having already its soution). In our context, the 
WKB parameter will be the inverse of the elementary volume. 

It is clear that with this method we can, at least in very simplified cases, make contact between the results of 
certain quantum gravity models concerning certain correlation functions and a language that is more close to the 
one of effective field theory. Clearly, what has been shown here is just what happens in the case of a minisuperspace 
model, in which all the complexity is reduced to the dynamics of a finite number of degrees of freedom. 

The hope is that the method can be generalized to the case of functional calculus (see, for instance, the detailed 
derivation of Einstein's equations by means of the Hamilton- J acobi method in |82j). and that it can be used, at 
least in principle, to develop a systematic effective field theory, for instance in the case of small perturbations from 
homogeneous configurations. To achieve this, however, we need more than what we have showed here. 



7 Corrections to the critical behavior 

In this section we are going to discuss briefiy the issue of the non-universal corrections to the critical behavior and 
their impact on the kind of observables that we have just mentioned. Again, as our case study, we will investigate 
the corrections to critical behavior of the form: 

W{g) « (g - gcV (l + 6i(g - Qc) + b2{g - 9cf + 0{{g - g,f)) (67) 

where 61,62 are the non-universal coefficients {i.e. model dependent) that control the first departure from the 
scaling behavior. 



7.1 Corrections to the partition function 

As in section [21 define 



In this case we get 

75c 



(.9 



1 dW 



+ (7 + ^i5c) + {bi - big, + 2625c) (5 - 5c) + 0{g - g^f (69) 



Notice that the presence of a non-vanishing 61 leads already to a correction to what would be the first term 
obtained expanding the results of section [2] beyond the lowest order. Working at the given approximation is it 
possible to find that: 

(5-5c)-- + ^ + ^ + 0(n-^), (70) 

where 

ai=75c; a2 = -ygch + bigc); ag = 75^7^ + 615^ + 6175^(8 - 6i.gc)] (71) 
As a consequence, quantities like f2(n) become, to the next to the next to the leading order, 

n{n) = - log gcU - 7 log n -I- ( - — -I- 7 log ai ) + — (72) 
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where ci is a calculable coefficient. Hence, taking the exponential, 

P{n) = exp (-an - ^) (73) 

In principle, one should insert this expression back into the machinery used to get the effective equations of 
motion. However, at large volume the effect of the logarithmic corrections due to the scaling (i.e. the prefactor 
vT') arc still larger than the non-universal corrections. On one side this is an interesting result, since it means that 
we could in principle detect as a first deviation from GR some operators that are associated purely to the critical 
behavior of the underlying model, while the corrections due to the granular nature of spacetime might become 
relevant at higher energies. On the other hand, this result might be shifting even farther away the access to the 
microphysics of the system, which is of course disappointing for the quest of the identification of the microscopic 
theory describing spacetime. 



7.2 Corrections to the wavefunction 

The corrections to the partition function are not enough to understand what is the real impact of the critical 
behavior on low energy effective theory. In the case of amplitudes computed with the given auxiliary fields, 
however, we need to include more into the game. Indeed, as anticipated in section 31 we should expect not only 
corrections due to the fact that the system is not exactly at criticality, but also the fact that the critical exponents 
might in turn depend on the position in the space {X,g). 
It is convenient to define 

e = 5-5c; V = >^~>^c, (74) 

where (A^, g^) is the critical point we are approaching^- In the neighborhood of this point we might Taylor-expand 
the various functions that are involved. For instance 

Ae(5)« A^ + rfe + 0(e)2, g,{X) ^ g^^ + crj (75) 
7(A) « 70 + a?? + O(ry'), S{g) ^ So + be + 0{e^) (76) 

whence 

Q(A, <?) = (A - A,(ff))^(s) (g - 5c(A))'^(^) (1 + m,7^ + n,e + 0{e\if,eii)) (77) 
might be approximated by 

Q{X,g) w (?7 - def°{e - cvi)''° (l + mir; + nie + 6elog(7/ - da) + aT]\og{e - crj) + 0{e'^ ,Tf- ,erii)) (78) 

The running of the critical exponents will lead to logarithmic terms which give rise to additional corrections to the 
wavefunction. Unfortunately we will not be able to work out an explicit analytical calculation which takes these 
terms into account. However, what is crucial to observe here is that the presence of such corrections will lead to an 
inevitable mixing between e and 77 in the definitions of the bulk volume nbuik and the boundary volumes riboundary- 

l^bulk — ^9 ! '^boundary — ^'^.-ia? V'^/ 

Q og Q OA 

and hence, when going for a representation of the probability amplitude in terms of these quantities, there will 
be contributions to the wavefunction which will take into account not only the continuum limit on the boundary, 
but also the continuum limit in the bulk, the corrections to scaling and the corrections to the critical exponents. 
Still, the logarithmic term coming from the prefactor will dominate over the corrections due to granularity, at long 
range. 



'^^Notice that, due to the presence of an auxihary group field theory, niight be difli^erent from the critical point for the partition 
function. 



17 



8 Dynamics of gravity and gravitational entropy 



One of the most intriguing aspects of GR that requires an explanation is the nature of the relationship between 
gravity and thermodynamics |83j . This relationship seems to be crucially related to the presence of horizons, in 
the broad sense of null surfaces, and the fact that general relativity can be rewritten in terms of laws describing 
the evolution of them. 

Clearly, this requires an explanation. Using the results of section [2] we now see how it is possible to provide 
a definition of entropy, by introducing a statistical distribution that describes the property of an ensemble of 
microscopic states, weighted by the thermodynamical potential n{n). 

On purely statistical grounds we can define a notion of entropy associated to this distribution function, 

/•oo 

5stat = -/ P{n)\og{P{n))dn (80) 
Jo 

To compute this expression we make use of standard manipulations: 

(81) 

rj=0 

We get: 

5stat = log S - (1 + 7 + 7 log(a)) + -f^f^^ = (82) 

^stat = logr(/3) - /3 log(r/3-i) - (/3 - 1) log(r/3-l) + $(/3) = -(2/3 + 1) log(y) + ^(/3) (83) 

where is a function of j3 alone. Remember that (3 controls the fluctuations: the larger is /?, the smaller are the 
fluctuations. Furthermore, /3 is controlled by the exponent 7 that describes the approach to criticality, and hence 
is a universal property of the system. This suggests to isolate all the dependence on the term 

The natural comparison that we can make is between the result ()83p and the partition function for euclidean 
quantum gravity with a cosmological constant, computed with the method of gravitational instantons in |85j . This 
comparison shows a basic discrepancy: we do not get the term scaling with the area of the cosmological horizon. 
This in a certain sense obvious, given that in this simple model we have not introduced any mean to control the 
contribution to each amplitude of the appropriate discretization of the Ricci scalar. The mismatch between the 
starting microscopic actions might be the first reason. A second reason could be that the critical behavior generating 
spacetime as we know it is of a different type. Unfortunately, knowing S'stat, which makes reference only on the 
average properties of the statistical distribution, does not allow us to go backwards and reconstruct the critical 
behavior of the partition function. While this term is per se intriguing (opening up a direct connection between 
macroscopic calculable properties of spacetime to microscopic ones associated to criticality), it is probably not the 
whole contribution to the gravitational entropy: it might be only the entropy associated with the uncertainty in 
the total volume. 

In fact, there is another statistical entropy that we can define, just observing the way in which the partition 
function is constructed in terms of triangulations F 

W^ = 5]g"(r)A(F) (84) 
r 

The function (7"'^'"^A(F) can be seen as a distribution function in the space of triangulations, and as such we can 
compute its statistical/Shannon entropy: 

Sw^~\ogW+^Y. ff"^'^'^(r) log (^"^^^^(r)) =-\ogW + n + J2 5"^^^^(r) log(A(F)) (85) 
r r 

This entropy would be much more compelling in light of the intuition that spacetime itself is a thermodynamical 
notion, result of the averaging of a lot of microscopic configurations. 



x^e '^log(a;''e '')dx 



dv Jo 



(x'' exp(-x))' 
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Unfortunately, we cannot yet devise, at the level of the microscopic action like the GFT one, a way to introduce 
a dummy parameter such that we can obtain this quantity by taking, for instance, a derivative of a generating 
function to get the desired result (as it happens, for instance, in the case of the evaluation of the entropy in the 
canonical ensemble as a suitable derivative with respect to the inverse temperature). What is clear from this 
expression, however, is that this entropy will be a certain function of 5 — ^c, or, equivalently, of n which is an 
extensive function of the average triangulation, in principle scaling with its volume. 

In fact, the main lesson that we can get out of the analysis of thermodynamics of gravity is that the thermody- 
namical properties are not localized on given regions, like in the case of black holes, but rather each infinitesimal 
region of spacetime, in its own, is a thermodynamical system 861. The case of black holes is just a situation in 
which it is particularly easy to uncover the relationship. 

However, functions like the entropies that we have defined so far depend on volumes, rather than surfaces, 
and tend to refer to global properties of the statistical system [i.e. of the classical spacetime) rather than to local 
properties as we should have to match the behavior of local Rindlcr horizons. 

There are at least three possibilities. First of all, the kind of formalism we have in mind is tailored for the 
case of Euclidean spaces, rather than Lorentzian spaces. Signature must be appropriately introduced, but also the 
dynamics probably needs to be revised, to get the appropriate scaling of entropies with areas and not with volumes. 
In this sense, the Lorentzian nature would be an essential building block. 

The second possibility is that horizon entropy is a sort of "red herring" : after all, it has been proven that horizon 
thermodynamics is not at all limited to the case of general relativity, and that many different geometrical theories of 
modified (but Lorentzian) gravity can be manipulated in such a way to reproduce some thermodynamical behavior, 
at the price of introducing non-equilibrium terms [87l HH [89l [90]. What is really important, in the pregeometric 
case, is that the amplitudes computed (and hence their semiclassical limit) match with some form of gravitational 
theory. 

The third case is in a sense mirror symmetric to the second one, and it is related to the fact that the kind of 
entropies that we are trying to define have their own dignity, but their definition is only formal, in the sense that 
they do not enter the balance of physical entropy encoded in the generalized second law (where all the forms of 
entropy, including the gravitational one, do enter the balance) . 

To settle this issue, one needs to introduce matter fields explicitly and hence define gravitational entropy in such 
a way that, together with matter entropy defined in the conventional way, it is possible to formulate consistently 
thermodynamics (with a generalized second law). This is a point of view dual to the one advocating the need to 
introduce matter fields for the very definition of the physical geometry. 

In this respect, there is another important (and related) issue involves the formulation of a low energy effective 
field theory for gravity: the implementation of local Lorentz invariance. Given the present constraints on Lorentz- 
violating effects in particle physics [91] . we are forced, if we do not want to introduce some additional microscopic 
symmetry yet to be discovered |921 1931 194j , Lorentz invariance at the fundamental level. Similar arguments have 
been shown to hold, in simplified cases of emergent gravity, for the equivalence principle |34| . However, one might 
still have some Lorentz breaking effects (spontaneous breaking as in [95], or anisotropic scaling features like in |56j). 

There is some interesting phenomenological implication that deserves to be studied. Within formalisms like 
GFT, the wavefunction and the correlation functions, even in the continuum limit, cannot be expected to implement 
locality at the level of spacetime, given that the starting point is nonlocal in an auxiliary space, and that spacetime 
emerges only in a statistical sense. Therefore, the effective action will be, generically, a complicated nonlocal object. 
The limit of local quantum field theory will have to emerge in some regime, with locality as implemented in the 
standard model being only a long range phenomenon. In this respect, it will be interesting to compare these models 
to other models of emergent locality discussed in other contexts of emergent /analogue gravity [23 , 24} [25 l [26 l [27 l 196 ) . 
as well as to develop a full phenomenological framework to test nonlocality effects in particle physics processes, by 
constructing effective field theories with nonlocality features. 

9 Concluding remarks and outlook 

In this paper we have tried to address the issue of the continuum semiclassical limit of class of models for quantum 
gravity like group field theory, basing the analysis on the conjecture that near criticality the behavior of the 
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connected part of the partition function is essentially that of a homogeneous function. Despite this restriction, that 
is due to criteria of simplicity rather than on a full fledged calculation which is still lacking, the general logic can 
be exported to more general case. 

The most significant step in the investigation is the introduction of suitable generating function for amplitudes 
or wavefunctions, in section|4l that allow us to use the critical behavior of the statistical system to infer the shape of 
the wavefunction of the universe in the sense of Hartle and Hawking, and hence, by techniques discussed in section[S] 
and [SI used to draw very primitive continuum effective field theories, at least in the case of minisuperspace models. 

Clearly, much more work needs to be done, but still the results we are collecting here seem to be promising, 
showing for instance that it is possible that the first corrections to the Einstein equations might tell us something 
about the continuum limit (its universality class, in particular), rather than the corrections due to the detailed (and 
largely unknown) microscopic structure of spacetime. This is of course much more interesting than the opposite 
situation, since it would allow us to discriminate between different microscopic models for spacetime just looking 
at their universality class. 

The procedure elucidates also the origin of the macroscopic coupling constants coming from the microscopic 
theory. While it is still premature to envisage a systematic identification of the various terms in a gravitational 
action in terms of curvature invariants, the general intuition is clear. They will be the outcome of the critical 
properties of the system used to define the continuum limit, with corrections coming from the granularity of 
spacetime (encoded in the non- universal terms). This is particularly significant since it sheds a different light on 
the problem of the renormalization of the gravitational coupling constants. For instance, the actual value of the 
cosmological constant will not be determined by any vacuum energy counting: there is no such a thing as an energy, 
in the pregeometric formalism we are using. The only thing that does matter is the way in which the continuum 
limit is approached, as well as the underlying dynamics. This is perfectly analogous to what has been discussed 
in [97] (see also [98l [99]). In light of this, the cosmological constant problem must be translated into a question 
regarding the critical behavior of the underlying theory. A pregeometric approach to the cosmological constant 
problem seems to be a rather intriguing one, deserving further investigations. 

While all these topics do have their relevance, a calculation that can elucidate many of the previous issues 
would be the determination of the statistical properties of De Sitter spacetime. De Sitter spacetime shares many 
thermodynamical properties with black hole spacetimes, given that cosmological horizons do define thermodynam- 
ical quantities associated to their dynamics [85j . However, the nature of these thermodynamical properties is much 
more interesting than the corresponding ones for black holes. Indeed, while one can imagine that the thermody- 
namical properties of black holes can be inferred by the properties of the degrees of freedom more or less localized 
within the black hole event horizon, the same cannot be said for De Sitter space. Indeed, cosmological horizons 
are observer dependent, given that the space is maximally symmetric. Hence, if the thermodynamical behavior of 
DeSitter space can be associated to some pregeometric degrees of freedom, they must be "equally distributed" in 
spacetime: in this sense, probing De Sitter spacetime we should be able to have more direct access to some of the 
properties of the fundamental degrees of freedom of spacetime [1031 1104] . 

As we have pointed out at the end of section [2l we do not have yet any positive result in this direction. In 
a scenario in which only equilateral simplices are glued (via a GFT-like formalism), it is impossible to include 
information about curvature, and hence to introduce even a very simple form of gravitational action. We need to 
go beyond this model. 

As we have shown, by introducing certain auxiliary fields we can work with amplitudes and wavefunctions, 
and thus infer from there the properties of the effective dynamics by a comparison with. However, we still lack a 
systematic way to control the boundary data in a local way: while it is easy to ask questions concerning homogenous 
spatial sections, it is unclear how to introduce concretely the inhomogeneities (and the parameters controlling them) 
associated to gravitons. 

Nonetheless, this result suggests that at least part of the logarithmic corrections to the entropy of de Sitter 
space can be originated by the statistical fiuctuations of the system near criticality, and hence lead to the conclusion 
that these corrections, if detected, can give direct access to the microscopic theory underlying classical spacetime, 
shedding some light onto the continuum limit. However, a reproduction of the De Sitter entropy calculation would 
not be enough. Indeed, what makes gravitational entropy relevant is that it enters the balance of entropy flows, 
and hence a big deal would be to check that the entropy as it would be defined is really counted in a generalized 
second law. 
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To conclude, let us examine briefly the other important limit that we have neglected so far, the semiclassical 
limit. From the results of sections [5] and [Slit is clear that, whenever we translate the results of the critical behavior 
of GFT in terms of a continuum field theory, it will still be quantum field theory 0, containing all the quantum 
corrections to the classical action. Therefore, we need more. Since we do not have free parameters to play with (an 
overall inverse temperature factor in front of the GFT action can be reabsorbed with a field redefinition) , we have 
to conjecture that the semiclassical limit will require some other large number to appear in the action. Again, for 
naturalness reasons, this number can be associated only to the presence of a large number of additional degrees of 
freedom, i.e. to a large number of matter fields. Obviously this conjecture must be carefully proven, but it would 
match some rather old ideas on the semiclassical limit of quantum gravity [100] . 

We want to conclude with one of the most puzzling points of this paper. In section [2] we have said that the 
models that we are playing with are leading to equilateral simplices, and that all the information about curvature is 
lost. Of course, this is a very unpleasant situation, at a first glance. However, the results of sections [5] and [5] clearly 
show that, despite this trivial short range dynamics, the long range dynamics, controlled by the properties of the 
critical point, might be nonetheless highly nontrivial, at the point of resembling the symmetry reduced dynamics 
of gravitational theories. 

This might be expected in view of the general argument of universality [1011 1102| : many different microscopic 
Hamiltonians, belonging to the same basin of attraction of a given critical point, will lead to indistinguishable 
critical behavior. This might be seen as an hint that perhaps, at this level of refinement, it is enough to construct 
a model with the correct universality class, rather than with the correct microscopic dynamics. 

This of course is still a difficult task, as the understanding of dynamical triangulations versus causal dynamical 
triangulations has told us. There, the inclusion at the microscopic level of all the relevant terms (the discretized 
Einstein-Hilbert action with cosmological constant) does not guarantee that the continuum limit will describe 
four dimensional spacetimes. Therefore, it might be well that the implementation of a microscopic GFT model 
might be complicated by the fact that the implementation of a simplicial path integral for gravity might be less 
important than the identification of the correct universality class, which is ultimately controlling the macroscopic 
limit. Therefore, the classification of GFT universality classes stands out as one of the most important problems 
that such a formalism has to face. 

On this possible mismatch between microscopic and microscopic level, the analysis presented in [431 about the 
effective actions around prescribed GFT vacua shows that, if these effective dynamics encodes the gravitons, their 
action might be quite different from the fundamental action with which we start building spacetime from nothing. 
This of course complicates the path towards a formulation of concrete recipes that can tell us how to build the 
microscopic theory once we know what is the macroscopic limit we want to approach. 

As mentioned in the introduction, this might be a coincidence. In fact, we can argue for the contrary. After all, 
the system we are playing with is still a system in which we are trying to realize a phase transition that will lead 
to some notion of continuum geometry. By taking the point of view of Landau, whatever is the microscopic theory, 
the macroscopic order parameter we should expect is a three geometry (in the case of the wavefunction) . The 
appearance of an effective geometrodynamics will be just a consequence of the fact that this order parameter will 
have to satisfy some equations of motion (Ginzburg-Landau equations), which are some form of geometrodynamics. 
Whether these equations are leading to general relativity or to some other theory is a more delicate matter, that 
will be carefully examined in future work. 
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Again, working in the Euclidean sector with the partition function does not change the result. 
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